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Abstract 

We analyse the structure of relativistic corrections in the inclusive produc- 
tion of polarized J/\E' from b-quark decay. The analysis is performed not only 
for the production channel in which the cc pair is a color sinlet, but also for 
the channels in which the cc is a color octet. We find that the correction 
in the color-singlet channel at the tree-level is completely determined by the 
decay constant of J/vl/, while in the color-octet channels the corrections are 
characterized by three matrix elements defined in NRQCD, whose numerical 
values are unknown. We discuss the impact of these corrections on the po- 
larized J/^-production, and the impact is so significant that the predictions 
based on the analysis for the considered process may be unreliable. Finally, 
we propose an integrated spin observable to measure the polarization of J/ty. 
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Quarkonium systems generally are thought as simpler than light hadrons and it may be 
easier to handle them in the framework of QCD. However, only recently we have been able to 
treat their decays and productions rigorously, based on a factorization with non-relativistic 
QCD(NRQCD) In the framework of the NRQCD factorization the effect of short distance 
is handled with perturbative QCD and the effect of long distance is parameterized with 
NRQCD matrix elements. The factorization is performed by utilizing the fact that a heavy 
quark inside a quarkonium moves with a small velocity v in the quarkonium rest frame 
and an expansion in v can be employed. In this factorization, inclusive productions of a 
quarkonium, e.g., like J/^f, can be imagined at the leading order of v as the following: The 
c- and c- quark are produced at a same space-time point, and then this pair is transformed 
into the J/^>. The production of the cc pair can be calculated with perturbative QCD, 
while the transformation is a nonperturbative process, which can be described by matrix 
elements defined in NRQCD. At higher order of v various effects are taken into account, 
e.g., relativistic effect, and the effect of the produced cc pair which does not have the same 
quantum numbers as those of J /if!. For J /iff the transition of a cc pair in a color-octet state 
happens at higher orders in the small velocity expansion, but in some cases, the produced cc 
pair is more likely in a color-octet state than in a color-singlet state, hence the production 
of J /ty through a color-octet cc pair is not negligible, and has an important contribution 
to the total production rate. Including this contribution one is able to fitting experimental 
results of the total production rate obtained in Tevatron [0-[§]. 

With the NRQCD factorization, not only the total production rate of a quakonium 
can be predicted, but also the polarization of the quarkonium, if it has a spin. Recent 
preliminary measurement || at Tevatron by CDF shows that the produced J /iff is polarized 
in the way which is unexpected from theoretical predictions. Several attempts to explain 
the discrepancy are made |||7|]. It should be kept in mind that these predictions are only 
based the theoretical analysison at the leading order in a s and at the leading order in v. 
For charmonia the velocity is not very small, typically v 2 ~ 0.3, and one-loop effects are 
generally substantial. Adding effects at higher orders of a s and of v, the predictions may 
be changed significantly. In this work we study the effects of higher orders in v on the 
polarization of J /iff produced in b-quark decay, to see how important the effects are. 

In the framework of heavy quark effective theory inclusive decays of a b-flavored hadron 
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can be thought approximately as inclusive decays of a free b-quark. The inclusive 6-quark 
decay into a J/f were studied at the leading order in v in ]§|||, in the one- loop correction 
from QCD for the unpolarized decay was studied. In these studies it is already shown that the 
color-octet cc pair plays an important role. In this work we will not only considered the color- 
singlet contribution but also the color-octet contribution. With the NRQCD factorization, 
the polarized decay width for the process 

b^J/^ + X (1) 

can be written as: 

r A (J/¥) = Y,C n (b - cc[n] + X)<0|O J /*[n]|0>, (2) 

n 

in which the coefficients C n (b — > cc[n] + X) describe the production of a cc pair in a state 
n, the matrix elements (0|O J /*[n]|0) characterize the transition of a cc pair in rz-state into 
the J/^f. The coefficients can be calculated with perturbative QCD because the production 
of a cc pair is a short- distance process, while the matrix elements represent nonperturbative 
effects and they are defined with operators in NRQCD. These matrix elements are scaled 
by the power of v with the rule of power counting in v ||1 1|| . The index A stands for the 



helicity of the J/^f, A = L is for the longitudinal polarization, A = T is for the transversal 
polarization. We will use the matching procedure proposed in ]]J to identify the operators 
at the next-to-leading order in v and to calculate the corresponding coefficients. 



To do the matching we need to consider the process 

b{p b ) -> c(pi) + c ( p 2 ) + X, (3) 

where the momenta are given in the brackets. With a Lorentz transformation we can boost 
the cc pair into its rest-frame and denote in the rest-frame the three- momentum of c and of 
c as q and — q respectively. The total decay width for the process in Eq.(3) can be written 

The matching condition reads: 

f(q) =Y,Cn(b - cc[n] + X)(0\O c ~ c [n}\0) . (5) 

n 
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In the above equation, the q-dependence in the right hand side is only contained in the 
matrix elements, in which the hadronic state is replaced by the partonic cc state. With 
Eq.(5) one can identify the matrix elements appearing in Eq.(2) and can determine the 
corresponding coefficients. 

The effective weak Hamiltonian for 6-quark decay is: 

= % E {V cb V: q [\cMcY{l - 7 5 )cg 7 ,(l - 75 )6 

V 2 q=s4 6 

+C s (fi)cT a r(l ~ 75)cgT a 7 M (l - 75)6]}- (6) 

We neglected the contributions of QCD penguin operators in H eS . T a (a = 1, • • • 8) is SU(3) 
color-matrix. The coefficients C\ and C 8 are related to the usual C± by 

d(/*) = 2C+(a*) - C_(/i), C 8 (jm) = C + (jj.) + C-(jj.). (7) 



With the one-loop evolution one can obtain 

Ci(m 6 ) 



0.18 (8) 



C 8 (m b ) 

This indicates that the cc pair is produced more likely in a color-octet state than in a color- 
singlet state. Hence the color-octet contribution is significant and should be included in the 
decay width. With the effective Hamiltonia it is straightforward to calculate T(q) and to 
determine the coefficients in Eq.(5). However, the power counting of v for matrix elements 
with partonic states is different than that for matrix elements with hadronic states. Before 
presenting our results, we briefly discuss the identification of operators in Eq.(2). In general 
the matrix elements in Eq.(2) take the form 

(0|O>] |0> = £<0|/C n |tf + X) (H + X\)C' n \0), (9) 
x 

where H denotes a quarkonium. K n and KJ n are operators defined in NRQCD. The general 
constraints can be obtained by considering symmetries of QCD or NRQCD. Because QCD 
respects charge conjugation- and parity symmetry, the operators must be transformed in the 
same way under C- and P transformation. The operators K, n and K' n can be classified with 
the weight j of £77(2) representations. Since the total angle momentum is conserved, the 
operators must have the same weight. By considering the approximated symmetry, spin- 
symmetry of NRQCD, further constraints can be obtained. With these constraints, one can 
identify the operators appearing in Eq.(2). 
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Now we give our results for T\(J/^). We write 

r A (j/vi/) = ri 1) + ri 8) , (10) 

where the index 1 or 8 stands for color singlet contributions or for color-octet contributions, 
respectively. The leading order in v for color singlet contributions is v°, while the leading 
order for color-octet contributions is v 4 . We expand the quantities: 



r (l) = r (l,,0) + r(1) „ 2) + 



v 



T W =T W) +T W) +(D ( v By (n) 

For color-singlet contributions at the leading order, there is only one matrix element repre- 
senting the transformation of a cc pair in 3 Si state into the J/^f. At this order we have: 
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where we used |V^| 2 + |V CS | 2 ~ 1. At order of v 2 , there is also one matrix element which 
represents the effect of the relative movement of the cc pair inside J/^> . The results reads: 

T * ~ ~^mV mb{mb ~ 4mJ m\ Gxi ^ b ] (13) 



with 



3 + 92y + 20% 2 + 576y 3 

G L {y) = 



6(1 + %) 
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2y(3 - Q8y - 304y 2 - 960y 3 ) 

GAy) = siTTW ' (14) 

The definition of the matrix elements (0|O 1 J/,I '( 3 Si)|0) and (0| J P 1 J/ *( 3 5' 1 ) |0> can be found 
in 0. These matrix elements can only be calculated nonperturbatively and they must be 
calculated with the same accuracy in order of v if one uses both to make predictions. 

For color octet contributions at the leading order in v, the cc pair can be in 1 S'o, 3 Si and 
3 Pi states, these states can be transformed into a J/$ through emission or absorption of 
soft gluons. Correspondingly, we have in the matching three types of matrix elements, these 
matrix elements can be reduced to three matrix elements with the symmetries mentioned 
before. They are: 
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= ^(AK(A)(0|0 8 J/,I '( 1 5o)|0) (15) 
^(0|xWV| J/V + X) (J/* + X|^W X |0) 

= ^(A)S*(A)(0|0 8 J/,I '( 3 5 1 )|0) (16) 
^(01x^(^7 x (-|D))V|J/* + X)(J/* + X|^TVx (-Id))'" X |0) 

if 

= ^(<M(A) • e*(A) - £ l (A)£*(A)(0|O8 / *( 3 Pi)|0) (17) 

In the above equations £j(A)(i = 1,2,3) denotes the polarization vector of «//\& in its rest- 
frame. The field ^(x) is the field in NRQCD for the c(c) quark. D is the spacial part of the 
covariant derivative . The definition of the operator O 8 7 ^*( 1 S , ), 0^( 3 Si) and 8 7/ '*( 3 Pi) 
can be found in JE]. The rotation invariance leads to Eq.(15), while spin symmetry is used 
for Eq.(16) and Eq.(17). The leading order of these matrices is at v 4 . It should be noted that 
Eq.(16) and Eq.(17) also hold at order of v 6 , because to keep the correct spin configuration 
the spin of the cc pair must be flipped twice through the interaction which violates the spin 
symmetry, this can only happen at order of v 8 . With the above identification we obtain the 
color-octet contribution at the leading order: 

^(8,v 4 ) ClG 2 F \V cb \ 2 2 2 2 

IV = 5 [m h — 4m„) 

L 288Tvm 3 b m c y b cJ 

■ {mKOlO^CSoM + m^Olof^OiO) + 8(0|O^*( 3 i\)|0>} 

p v * } = ClGl\V cb \* 2 _ 22 

• Ll(0\OrCSo)\0) +4ml(0\Or( 3 S 1 )\0) + ^^<0|O 8 J/ *( 3 ^)|0>) . (18) 

mi 



At order of v 6 for the color-octet contributions various matrix elements appear at the 
first look. For example, the matrix element 

£(0| x W(-^) 2 D ( D*V| J/* + X)(J/y + X|^rV X |0> + h.c. (19) 
x 1 

To take a close look at the matrix element we consider 

B ijlk = ^( | x t T v i (-^) 2 D ( D^Vl^/^ + X)(J/V + X\^T a a 3 X \0) + h.c. (20) 
x 2 
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At order of v & the spin symmetry still can be used for the matrix element, it leads: 

& lk = e i {\)e* j (\)-A lk . (21) 

By rotation invariance, A ik is proportional to 8^- Because the tensor D D is symmetric 
and trace-less, we conclude that B' ljlk is zero at this order, hence also the matrix element in 
Eq.(19). This can also be understood as the following: Because the spin symmetry holds, 
the total orbit angle moment is conserved. The X-state in the bra in Eq.(20) is a / = 
state because of the conservation, while the X-state in the ket is a I = 2 state. Therefore, 
the matrix element is excluded by the definition of the sum over X state. By considering all 
constraints, only three types of matrix elements remain, and they can be further reduced to 
three matrix elements: 

EM^n-^VlJ/* + X)(J/* + X\^T a X \0) + h.c. 
x 1 

= ^(AK(A)(0|P 8 J/ *( 1 5 )|0) (22) 

]r(0|x f TV(-^D)V| J/* + X)(J/V + Xl^TVxlO) + h.c. 
x 1 

= ^(^(AXOlPg^^OlO) (23) 
£<0| X tT> x (-1d))*(-^D)V| J/V + X)(J/V + X^TV x (-^% X |0> + h.c. 

= ^(<Jij-e(A) • e*(A) - £;(A)£*(A)(0|P 8 J/ *( 3 Pi)|0). (24) 

The matrix elements (0|P 8 J/,3/ (%)|0), (0|P 8 J/ *( 3 Si)|0) and (0|P 8 J/ *( 3 Pi)|0) can be obtained 
by summing over the helicity and by the contraction over the indices in the above equations. 
With these matrix elements the results for the color-octet contribution at order of v 6 are: 

r (s,v 6 ) C%G F \V cb \ 2 2 2 
r A = 2887rm a m b (m b -Am c ) 

■ {F{^) ■ <0|P 8 J/ *(%)|0> + G A (4) ' (0|P 8 J/ *( 3 ^)|0) 
m b m b 

+^(^|)-i(0|P 8 J/ *( 3 P 1 )|0)} (25) 



ml m 2 



with 



H L (y) 



7 + 10% + UAy 2 + 320y 3 
6(1 + Ayf ' 
%(1 + My + 240j/ 2 + 704y 3 ) 
6(1 + 4y) 2 ' 
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= _7 + 84 y +14V + 70V 
vy; 6(1 + 4|/) 2 v ; 

The function G\(y) is given in Eq.(14). The above results indicate that the correction in the 

color-octet contributions is to take the effect of the relative movement of the cc pair inside 

the J/^S into account. This is similar as the case with the color-singlet contribution. Hence, 

the total corrections are just relativistic corrections characterized by four matrix elements. 

With the results given above one may predict the polarization of the produced J/^> and 
the total decay width T = 2IY + if one knows the numerical values of the eight matrix 
elements. Among the four matrix elements at the leading order in v the best known matrix 
element is (0|Of ( 3 5i)|0), which is calculated with potential models and with lattice QCD, 
while the other three are extracted from experimental data. The value of {0\Og^ ( 3 Si)\0) 
is rather well determined by direct J/ty production at large transverse momentum in pp 
collisions [0-fl. However the uncertainty for this determination is large and it can be at the 
level of 100%. The other two are not well determined, only certain combination of them is 



know with a large uncertainty. In ]TJ| a re-analysis of experimental data from Tevatron and 
from Hera is performed, in which some higher-orders effects due to multiple-gluon initial- 
states radiation are considered. It is shown that the values of these matrix elements can 
be changed substantially by including these effects. Based on the leading order results an 
analysis for predictions of the polarized J/\I/ is given [[|. The four matrix elements at the 
next-to- leading order in v are completely unknown. With the power-counting in v, we only 
know that they are suppressed by v 2 relatively to the corresponding matrix elements at the 
leading order in v. All of these prevents us from a detailed prediction for the polarization 
and for the total decay width. Nevertheless, one can still see the impact of these corrections. 
For the color-singlet contributions, if one neglects the one-loop QCD correction and uses 
the vacuum saturation, then with H e ^ in Eq.(6) one has only one constant representing the 
nonperturbative effect in the process in Eq.(l): 

(J/V\crc\0) = -ifj/vMj/vie^X))*, (27) 

where fj/^ is the leptonic decay constant of J/^f and is related to the leptonic decay width: 

r(j/tf -> £ + r) = cm f 2 j^. (28) 

In the above equation, the only approximation is to neglect effects of higher orders in a em and 
lepton masses, relativistic corrections are automatically included. In NRQCD the vacuum 
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saturation brings uncertainty at order of v 4 . Therefore the color-singlet contributions can 
be written as: 

rilril It/ 12 

,(1) _ W^Fl v cb\ ( 2 T,/r2 ^2*^2 r2 , /n( 2\ , fi^( n A\ 



if = C ;^t' 2 K - ^/*) 2 "W/» + 0(a s 2 ) + 0(A (29) 



where the relativistic correction calculated before is contained in the J/\l/-mass Mj/^ and 
in the decay constant. Comparing the results in Eq.(29) with that in Eq.(12) we may see 
how large the relativistic correction is. For this we define: 

F (M°) 

*P = TOT' (3°) 

where T^' v ■* is given in Eq.(12) and the results in Eq.(29) are used for rf . To obtain 
numerical values of we take = 4.7GeV and m c = 1.5GeV. fj/^ is determined 

by the leptonic decay to be 405MeV. For the matrix element (0|O^*( 3 Si)|0), its value is 
determined based on a potential model in |HJ to be 1.16(GeV) 3 , which is in agreement with 
a calculation of lattice QCD in |Tj|. With these values we obtain: 

« R ( t ] w 1.6. (31) 

This indicates that the correction is negative and very large. Another way to see the impact 
of the correction is to compare the coefficients in the front of the matrix element. We take 
the above numerical values for quark masses and obtain: 

r (D = CfGlW . U mo J^ Sim _ 52.4^<0|P 1 J/ *( 3 S 1 )|0>) + Q(v% 
43271 mi 



r g) = 9l9^L . J 9 .90(0|O 1 J/ *( 3 S 1 )|0) - 11.4-^(01^(^)10) \ + 0(v% (32) 
4327T | mi J 

With these results the relativistic correction can be at the level of 60% for 1^ and at the 
level of 35% for 1^ if one takes 



— 2 (0\P 1 J/ "( 3 S 1 )\0)^v 2 (0\O( / "'( 3 S 1 )\0), v 2 ^0.3. (33) 
mi 

Similarly we obtain for the color-octet contributions: 

r ( 8) = C iG 2 M 2 . {24.4(010^(^0)10) - 68.3 ■ -^<0|P 8 J/ *(%)|0) 



c 
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+24.4(0|O 8 J/ *( 3 S 1 )|0> - 52.4^(0|P 8 7/ *( 3 S 1 )|0) 

m c 

+ 19.8^(010^(^)10) - 35.8^(01^(^)10)} + 0{v% 
r (8) = c l G Vy^ . {24.4(0|O^( 1 i S )|0> - 68.3 • i^P/^ 1 ^) |0> 



2887T x . » v »/. / m 



c 



+9.9(0|O 8 J/,3/ ( 3 5 1 )|0) - 11.4^(0|P 8 J/ *( 3 ^)|0) 



+34.3^(0|O 8 J/ *( 3 P 1 )|0) - 86.5^<0|P 8 J/ *( 3 A)|0>} + 0(v 8 ). (34) 



If one assumes that there are similar relations among the color-octet matrix elements like 
that in Eq.(33), then the correction is very large. For example, the correction for the 1 S'o 
production channel can be at the level of 80%, and the correction for the 3 Px production 
channel can be at the level of 70%. With the above discussions one may conclude that the 
predictions based on the leading- and next-to-leading order in the small velocity expansion 
for the process are unreliable. 

The last subject of this work is to propose an integrated observable to measure the spin 
of the produced J/ty. Usually, one looks at the leptonic decay of J/ty to measure the spin. 
At a fixed momentum P of J/\l/ one measures the distribution: 

-(J/* -> tC) oc 1 + acos 2 #, (35) 

a cos 

where 9 is the angle between P and k, k is the momentum of the lepton in the J/^-rest 
frame. The parameter a is predicted as: 

a = Tv+TV (36) 

This may have a disadvantage that it will be hard to determine the distribution if the number 
of events with a fixed P is small, hence the parameter a. We propose to use an integrated 
spin observable to overcome this disadvantage. For this purpose we define the density matrix 
Rij of the produce J/^fr. For arbitrary polarization the decay width can be written: 

r(6 -> J/tf + X) = -i- / dne i R lj (P)e 3 , (37) 

47T J 

where Q is the solid angle of P, P denotes the direction of P. Similarly we can define the 
density matrix for the leptonic decay. Any observable O can be then predicted by 

^^/S/fr '^-^ (38) 
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where N is a normalization factor so that (1) = I. A simple calculation leads: 

Pij (k) = ^ - ^(kkj - ^Sij), (39) 

where py is normalized. Rij can be written: 

Rij(P) = a5ij + ibe ijk P k + cPiPj, 

a = r T , c = T L - r T . (40) 

The anti-symmetric part exits due to that parity is violated, and it is irrelevant here because 
the parity is conserved in the leptonic decay. With the tensor structure one can construct 
the integrated spin observable Op, and predict its value with Eq.(38): 

P = (P.kf-i, (O f ) = ggz__l£) . (41) 

If (Op) =0 the J/ty is unpolarized. If one know the values of the matrix elements discussed 
before, one can obtain the value for (Op) to compare with the measured in experiment. One 
may also obtain the invariance (Op) to determine the statistical error of (0 P ). 

We summarize our work: We have analyzed the relativistic corrections for the polarized 
J/^-production in 6-quark decay. We have calculated the perturbative coefficients and 
identified the matrix elements at the next-to- leading order in v. We find that the corrections 
can be very large in the color-singlet production channel and as well as in the color-octet 
production channels. For the color-singlet production channel the correction is determined 
by the leptonic decay constant of J/ty. These corrections are so large that the predictions 
for the considered process may be unreliable, if one only keeps several leading terms in the 
small velocity expansion. An integrated spin observable is proposed to measure the J/^f 
polarization. If detailed information of the eight matrix elements is known, numerical values 
of polarized decay widths and the observable may be obtained. 
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